Astronomy & Astrophysics manuscript no. ArtGRBdefinitif 


© ESO 2009 


January 22, 2009 





Luminosity of a quark star undergoing torsional oscillations 

and the problem of 7-ray bursts 

J. Heyvaerts^, S. Bonazzola\ M. Bejger'-^, P. Haensel^ 
' LUTh, UMR 8102 du CNRS, PI. Jules Janssen, 92195 Meudon, France 

- N. Copernicus Astronomical Center, Polish Academy of Sciences, Bartycka 18, PL-00-716 Warszawa, Poland 
^ Observatoire Astronomique de Strasbourg, 1 1 rue de I'Universite, 67000 Strasbourg, France 

Received 16 November 2007 / Accepted 5 December 2008 

ABSTRACT 

Aims. We discuss whether the winding-up of the magnetic field by differential rotation in a new-born quark star can produce a 
sufficiently-high, energy, emission rate of sufficiently long duration to explain long gamma-ray bursts. 

Methods. In the context of magnetohydrodynamics, we study the torsional oscillations and energy extraction from a new-bom, hot, 
differentially-rotating quark star. 

Results. The new-born compact star is a rapid rotator that produces a relativistic, leptonic wind. The star's torsional oscillation mod- 
ulates this wind emission considerably when it is odd and of sufficient amplitude, which is relatively easy to reach. Odd oscillations 
may occur just after the formation of a quark star. Other asymmetries can cause similar effects. The buoyancy of wound-up magnetic 
fields is inhibited, or its effects are limited, by a variety of different mechanisms. Direct electromagnetic emission by the torsional 
oscillation in either an outside vacuum or the leptonic wind surrounding the compact object is found to be insignificant. In contrast, 
the twist given to the outer magnetic field by an odd torsional oscillation is generally sufficient to open the star's magnetosphere. The 
Poynting emission of the star in its leptonic environment is then radiated from all of its surface and is enhanced considerably during 
these open episodes, tapping at the bulk rotational energy of the star. This results in intense energy shedding in the first tens of minutes 
after the collapse of magnetized quark stars with an initial poloidal field of order of lO''' gauss, sufficient to explain long gamma-ray 
bursts. 

Key words, gamma-ray bursts - quark stars - magnetic field - plasma - stars 



; 1. Introduction and motivation 

Understanding the physical nature of y-ray bursts (GRBs) in a 
way that is consistent with observations of their entire evolution 
remains a challenging mystery. A vast literature on the subject 
■ exists and a large number of models have b een proposed to ex- 
' plain this phenomenon (fo r a review see e.g.. lZhang & MeszarosI 
! 120041 and lM6szarosl200 6'). In this paper, we limit ourselves to the 
' case of long GRBs, of duration of between about 10 s and 1000 
' s. There are several basic facts to be explained. First, the re- 
lease of about 10"*^ - 10^' ergs in y-rays, of mean power higher 
than 10"^^ ergs s"'. The violent energy outflow, which eventu- 
ally transforms into a GRB, must originate in a compact volume 
of linear size ~ 10^ cm, because of the observed millisecond 
variability of the GRBs. To achieve the observed bulk Lorentz 
factor r = 100-1000, an energy outflow of 10^' ergs should 
have a rest-mass load of only 10^^-10"* M©. In other words, 
the baryon wind associated with long GRBs is only 10 '^-lO"* 
M0 s Therefore, within the inner engine of GRBs the sepa- 
ration of light from the matter is realized, producing the most 
luminous electromagnetic explosions in the Universe. 

Quark stars are hypothetical stars that consist of deconfined 
quarks (the s tructure o f these stars was studied in detail by 
IHaensel et alJil986 and lAlcock et al.iri986a) . They are presum- 
ably born in special supernovae, such as SNIc, from the col- 
lapse of very massive Wolf Rayet stars dPaczviiski & Haensell 
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I2OO5I) . The quark star itself forms in a second collapse a few 
minutes after the new-born proto-neutron star simultaneously 
deleptonizes and spins up. Alternatively, a quark star could re- 
sult from the collaps e of accreting neutron stars in X-ray binaries 
(ICheng & Dailll996h . The collapse of stars of an initial mass less 
than 30 M© is expected to resul t in th e formation of a compact 
object rather than a black hole dFrver & Kalogera 2001). Some 
authors claim that this limit is in fact 50 (IGaensler et^ 
120051) . which implies that a large fraction of the progenitors 
of SNIc's should eventually evolve into compact stars. Among 
those, a possibly non-negligible fraction could be quark stars 
(Paczyiiski & Haensel 2005). 

Bare quark stars differ from the normal, nucleonic, neutron 
stars in that the surface of a bare quark star (a strange star) plays 
the role of a membrane from which only leptons and photons 
can escape. This property was noted already in the early papers 
about quark stars d Alcock et alJl986al:lHaensel et aljl991 ). The 
quark star surface then eff'ectively separates baryonic from both 
leptonic matter and radiation. The exterior of a quark star is free 
of baryons, since the latter would be accreted onto the star and 
converted into deconfined quarks, resulting in a release of en- 
ergy. The close environment of a newly-born quark star is there- 
fore expected to be baryon-free. 

If quark stars do indeed exist, they would be prime candi- 
dates for emitting relativistic winds with small baryonic pol- 
lution. The importance of small baryonic pollution in the 
context of GRB fire ball models was already emphasized by 
IHaensel et al.l dl99ll) . when discussing the collision of quark 
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stars as an inner engine of short gamma-ray bursts. Models of 
the emission of gamma radiation in a GRB (De Ruiulal ll987l: 
[Zhang & Meszaros 2004; Dar & De Riijula 2004) require that 
the central engine emits a bulk relativi stic outflow with a Lorent z 
factor r ranging from 100 to 1000 (Meszarosll20"06l:lDaijl2006l) . 
Such a high Lorentz factor can only be reache d if the bary o n con- 
tent of the outflow is very low. For example Paczy risldl (1 19901) 
has shown that a radiation-driven wind can reach a Lorentz 
factor of 100 only if the luminosity injected in the wind ex- 
ceeds b y a factor 10^ the rest mass energy blown away per 
second. iBucciantini et al.l (l2006l) indicat ed (using th e results of 
ID calculation of relativistic winds by Michel 1969) that low 
baryonic pollution is necessary to obtain h igh Lorentz factors 
in a centrifugally-driven magnetized wind. iDessartetal.] 120071 
reached the same conclusion. For this reason, the quark star 
model is preferable because the low baryonic pollution of a 
quark star's environment ensures that energy can be deposited 
cleanly outside the star in the form of accelerated electron- 
positron pairs and y ray radiation without unnecessarily accel- 
erating baryons. This is why we strongly favour a quark star 
model and consider that the central engine of a long GRB may 
be such a star. The fact that quark stars could be the source of 
GRBs was first suggested bv lAlcock et al.l(ll986bl) . Our calcula- 
tions are however not specific to a quark star, except in Sects. [Til 
and l3.2l They would also apply to a strongly magnetized neutron 
star. 

The way in which the relativistic wind energy from the quark 
star is eventually released as gamma ray radiation is model- 
dependent and could be due to a range of environments at a dis- 
tance far larger than the light-cylinder radius from the star. Low 
baryonic pollution is needed only within a few light-cylinder 
radii from the central engine, a few thousand km, where the 
magnetized relativistic wind is expected to be accelerated to the 
required high Lorentz factors. Since the collapse from proto- 
neutron star to quark star is slightly delayed, no thick envelope 
is expected to be present in the rather small region where the 
leptonic wind is accelerated. 

Newly formed quark stars could be at the origin of 
GRBs because of the sudden transformation of hadronic 
matter into deconfined quark-matter when a neutron star 
or a proto-neutron star collapses to a quark star. Diff'erent 
ways in which such a collapse could be triggered have 
been described by a number of authors (Cheng & Dai 1996', 
1998a; Dai & Lu 1998a; Cheng & Dai 1998b; Bombaci & Datta 
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2002t lOuved & Sannin o"2002; Berezhia ni et al., ,2002> ,20031 
Drago et al.l l2004allbl; iP aczvhski & Haensell |200'1 iDrago etall 
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Another promising subclass of models for a GRB central en- 
gine is based on the existence of a rapidly (millisecond period) 
and differentially rotating compact star endowed with a strong 
(10'^ to 10'^ gauss) surface magnetic field, or spontaneously de- 
veloping such a field from differential rotation. This millisec- 
ond mag netar model of the GRB central engine was first pro- 
posed by lUsovl (Il992h . who suggested that the GRB energy re- 
lease derived its origin from the pulsar activity of a millisecond- 
period compact object with a dipole field of 10'^ gauss. The 
idea of a millisecond magnet ar has been revisited and dis- 
cussed by m any authors ( Thompson' ' 1 994^; ^Z hang & MeszarosI 
1200 IL I2OO2I) . Following Kluzniak & Rudermaiill998t), a num- 
ber of authors proposed differential rotation as a mechanism 
for strengthening the toroidal m agnetic fiel d in the interior of 
a newly-born neutron star (Ru derman et al.ll2000i) or in an ac- 
creting neutron star that develops difi'erential rotation as a re- 



sult of an r-mode instability (ISpruitlll999l) . Amplified magnetic 
fields, of the order of 10'^ gauss, would then be brought to the 
star's surface by buoyancy forces where the e nergy would be 
emitt ed in the form of bursts and Poynting flux. Ruderman et al.l 
(I2OOO) suggested that the emergence of strong fields would gen- 
erate episo dic, pulsar activit y from the open magnetosphere of 
the object. iDai & Lul (|1998b) considered this model for quark 
stars. However, it would be interesting to identify mechanisms 
that could generate the GRB phenomenon in objects with fields 
weaker than 10'^ gauss. In this paper, we propose such a mech- 
anism. 

A newly-born, compact star is entirely fluid. It supports dif- 
ferential rotation and internal oscillations of large amplitude that 
act on its internal magnetic field. In particular, differential rota- 
tion in the star would create a toroidal field from the poloidal 
component. This process is often referred to as an co process, or 
winding-up of the magnetic field. The wound-up field and the 
differential rotation velocity constitute an energy reservoir into 
which electromagnetic emission could tap. More generally, the 
star's internal motions may have some effect on the Poynting 
power radiated. The main problem of the quark-magnetar model 
is to explain how the existence of the wound-up magnetic field 
could affect the star's emission, either by direct extraction of the 
associated energy or by any other means. 

Our aim in this paper is to discuss the efficiency of the vari- 
ous ways in which a differentially-rotating, magnetized, compact 
star could shed energy in its environment. 

We restrict our attention to the case of aligned rotators in 
which the magnetic dipole axis is parallel to the rotation axis. 
The magnetic-field distribution in the star, although it may be 
locally structured, is regarded as smoothly distributed on a larger 
scale. 

We first show that field winding-up by differential rotation 
is not a steady, but oscillating process (Sect. |2]i. When much 
of the available kinetic energy of the differential rotation has 
been transformed into toroidal magnetic energy, the magnetic 
tension forces react back and reverse the motion. By means of 
this mechanism, any initial differential rotation develops into a 
torsional, standing, Alfven wave in the star. In the absence of 
losses, the star would oscillate, in the rest frame accompanying 
the average rotation, l ike a torsion pendulum (see for instance, 
iBonazzola et al.ll2007b . In Sect. |2] we determine the amplitude 
of the wave and in Sect. I3.1l we discuss its damping in a quark 
star as a result of second viscosity. This damping is found to be 
small. 

We then discuss various mechanisms by which the energy of 
the torsional oscillation could emerge from the star. Buoyancy 
is a possibility. It may however be quenched if, while the 
buoyant matter moves, some of the weak reactions between 
quarks remain frozen, which occurs if the matter is in a colour- 
superconducting state with a large enough gap, and the strange 
quark is sufficiently massive (Sect. 13). Alternatively, buoyancy 
may be inhibited by magnetic stratification or, if it develops, it 
could only redistribute flux in the star if the latter is magnetized 
in bulk. 

Since the internal magnetic field is time-variable, it could 
conceivably act as an antenna and radiate a large-amplitude, 
electromagnetic wave in an external vacuum. We calculate this 
radiation in (Sect. 14. Il l and find that the emitted power is insuf- 
ficient to produce a GRB. Radiation in a leptonic medium sur- 
rounding the star is shown to be equally inefficient (Sect.|4|2]l. 

We next consider the modulation by differential rotation of 
the DC Poynting emission of the fast, aligned magnetic rotator 
(Sect. 15). We find that an even oscillation, in which the south- 
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ern hemisphere oscillates in phase with respect to the north- 
ern hemisphere and the magnetic structure has a dipolar-type 
of symmetry causes a negligible modulation of the energy out- 
put. However, an odd torsional oscillation, in which the south- 
ern hemisphere oscillates in phase opposition with respect to the 
northern hemisphere, easily causes the star's magnetosphere to 
be blown open in a time-dependent way. An even oscillation act- 
ing on a magnetospheric structure that would not be strictly sym- 
metrical with respect to the equator has a similar effect. These 
openings of the magnetosphere cause a modulation of the power 
emitted in the relativistic wind blown by the fast rotator that is 
large enough to meet the requirements of energy and time scale 
necessary to explain the GRB phenomenon, even for moderately 
magnetized stars, with a field of order of a few 10'"* gauss (Sects. 
I5.3ti5.4l ). A collapse that is strictly symmetrical with respect to 
the equator would not excite odd oscillations. However, the ex- 
istence of a kick received by neutron stars at their birth indicates 
that a supernova collapse is in reality not strictly symmetrical. 
We show that even a weak, odd oscillation is sufficient to open 
the star's magnetosphere during several tens of minutes after the 
collapse. 

We use a Gaussian CGS system of units throughout the paper 
and spherical polar coordinates based on the rotation axis, r, 0, 
and (p, where 6 is the colatitude. The corresponding unit vectors 
are e^, eg, and e^. 



2. Torsional oscillation in the star 

Differential rotation necessarily induces, in an highly magne- 
tized and conducting star, a torsional oscillation. Such Alfvenic 
oscillations in fluid, magnetized, com pact stars were studied by, 
e. g., Bastrukov & Podgainy ( 1996), Rincon & Rieutord (2003), 
and Bo nazzola et al.l (12007). We assume that inside the star the 
poloidal part of the magnetic field is time-independent. This is 
reasonable because the strange matter is but weakly compress- 
ible. We also assume that perfect MHD is valid. 

2.1. Magnetic diffusion time scale 

Perfect MHD is a good approximation when the magnetic difiii- 
sion time tb is longer than the timescale of the considered phe- 
nomenon. The time tb depends on the electrical conductivity 
cr^ and the gradient lengthscale of the field which we assume to 
equal the star's radius R, such that tb ~ A-no-gR^/c^. The electric 
conductivity is the sum of the electronic and quark conductivity, 
(e) and cTeiq). The electron fraction in quark matter depends 
considerably on the physical conditions. The quark conductiv- 
ity wa s calculated for normal quarks by iHeiselberg & Pethickl 
(Il993h . Dynamical screening of transverse interactions by the 
Landau damping of the exchanged gluons is important in deter- 
mining quark mobility. The result can be expressed, for normal, 
massless quarks and a strong coupling constant supposedly equal 
to 0.1, as: 



cTAq)^5.6\0'^T-'l^ (fihlnofl^i 



(1) 



where ni, is the baryonic number density, no the nuclear den- 
sity (0.17 fm-^^), and Tio = r/10'°°K. If quarks are colour- 
superconducting in a colour-flavour locked (CFL) state, charge 
neutrality of the quark component alone is enforced because 
pairing is maximized when nj - «„ - iis. This happens when 
the mass m, of the strange quark is insufficiently large com- 
pared to the gap A. In the absence of electr ons, the CFL colour- 
superconductor is an electric insulator jAlford et al.l l2007l) . 



Electron suppression results in a difference between the chem- 
ical potentials of s and d quarks, which must remain limited. 
Electrons are suppressed tRaiagopal & Wilczek,2001,) only if: 



4/7 



< A, 



(2) 



where /i is a mean quark chemical potential. When the mass of 
the strange quark is so large that the inequality (|2]) is not sat- 
isfied, electrons remain present in the colour-superconducting 
quark matter and this matter is then a conductor with a conduc- 
tivity which may be far larger than Eq.([T]). If quarks are in a two- 
flavour colour-superconductor (2SC) state, only d and u quarks 
of two colours are paired, and electrons are present in quark mat- 
ter Then cri,{q) differs only by numerical factors from Eq.([T]| and 
crp(e) does not vanish. The magnetic diffusion timescale tb as- 
sociated with the scale length R and the conductivity ([1]) is about 
8 X 10'" seconds for iij, - «o and Tiq - 1. This is a lower 
bound to the true magnetic diffusion timescale, which is longer 
than this value when the electronic conductivity does not van- 
ish. This time is long enough for perfect MHD to be valid on the 
timescale of a burst. The exception to this general conclusion is 
when quarks are in a colour CFL-superconductor state in most of 
the star, and the s quark is so light that inequality (|2]) holds true, 
in which case quark matter is an insulator The star's interior is 
then electrically inactive and none of the phenomena described 
in Sects. ( 12.21) and (15.3b occurs. In this paper, we assume that the 
star's interior is a good electrical conductor, which implies that 
for the electrical conductivities implied by Eq. ([1]) perfect MHD 
holds true. 



2.2. Period and amplitude of the torsional oscillation 

In perfect MHD, the evolution equations for the velocity V and 
for the magnetic field B inside the star can be written, in a 
Galilean rest frame, as: 

^+iV-V)V-—(VAB)AB + V(£+U)^0, (3) 
of A-np 



^-VA(yAB) = 0. 

ot 



(4) 



where p, e, and U are the baryonic mass density, the enthalpy 
(the fluid is supposed to be barotropic), and the gravitational po- 
tential, respectively. The equations for the toroidal components 
of V and B can be written in the form : 



dt 4np r sin 
dB 



l—Bp-V(rsm9BA 



— - = rsin^Bp ■ V, 
ot \r sin 6 



(5) 
(6) 



We note that for rigid rotation, i.e. for Vsir, &) - Q,r sin0 with 
constant Q,, the r. h. s. of equation (|6]l vanishes, so that no 
winding-up of the field in this case occurs. Winding-up results 
from differential rotation with depth or latitude, or both. The 
Proudman-Taylor theorem does not apply when non-potential 
forces, such as the magnetic-tension force, are exerted onto the 
fluid and the motions are non-stationary. It is however possible 
that when non-potential forces are small compared to pressure 
forces and the timescale of the internal motions is long com- 
pared to the star rotation period, these motions organize them- 
selves such that becomes a function of the axial distance 
r sin 6 only. In that case, the field winding is only possible when 
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the poloidal magnetic field has a component perpendicular to the 
rotation axis. The buildup of the toroidal field causes a reac- 
tion magnetic-tension force (the terms on the right hand side of 
Eq. to grow. A torsional oscillation then develops. The sys- 
tem of Eqs. (|5]l and (|6]l for the azimuthal velocity and field inside 
the star has b een solved numeric ally with appropriate boundary 
conditions ( Bonazzola et al.ll2007 ). 

We define R to be the star's radius and Ru) - R/ilO km). 
The period Pj of the torsional oscillation is esssentially that of 
an Alfven wave with a node at both poles, that is, with a wave- 
length 2nR in a poloidal field Bp - Bpu lO''* gauss. We adopt 
lO''* gauss as a reference value since fields of order of between 
5 10'^ and lO'^ gauss are commonplace in isolated neutron stars 
(lHaberlll2007h and fields of several lO''* gauss are reported to be 
typi cal of anomalous X-ray pulsars and soft gamma-ray bursters 
dZLolkowski 2002). We assume that the field is initially rooted 
deep insi de the star. This view is sup ported by simulations of 
collapse (lObergaulinger et al] l2006al) . which indicate that the 
magnetic field after collapse is concentrated in the inner core. 
The density of the medium inside the star is of the order of the 
mean star density p» - M/V, M being the mass of the star and 
V = 47tR^/3 its volume. The period Pj of the torsional oscilla- 
tion is that of an Alfven wave of wavelength InR in a medium 
of density p.: 



Pt = 



BlR 



= 4.9s 



(7) 



This result may also be obtained from a linearization of Eqs. 
(|6|. The magnetic field, instead of pervading all of the star, could 
conceivably be present only in some superficial layer where the 
density is p < p», although the aforementioned simulations do 
not support this. In this case, Pj would be smaller than given by 
the expression in Eq.© by a factor ^Jp/pt■ The surface density 
of a quark star is 4S/c^, where S is the bag constant. For S - 
40 MeV fm"^, this surface density is 2.8 x lO^'* g cm""*, which is 
about 0.6 times the average density p.. For comparable poloidal 
fields, the period of a superficial torsional wave is only slightly 
less than the period given by Eq.© of a global torsional wave. 

The wave amplitude is set by the amount of energy Wd ini- 
tially stored in the collapse as kinetic energy of the differential 
rotation. We assume Wd to be a fraction ao of the order of a 
fe w percent of the star's total rotational energy, W*. Simulations 
bv lObergaulinger et alj (l200 6a b) indicated that ap does not ex- 
ceed 1 0% after the collapse to a neutron star. iBurrows et alj 
( l2007h found that the value of qd is less constrained. We adopt 
ao = 0.1 as a representative upper limit. We define /» to be 
the moment of inertia of the star and 745 = /^/(lO'*^ g cm^), P» 
to be the star 's average rotation period, an d Pt(ms) its value in 
milliseconds. lObergaulinger et alj (l2006albt) found that when the 
collapsed core reaches quasi-equilibrium, P^ most often ranges 
in value between 5 and 40 milliseconds. iBurrows et al.l (l2007h 
found that 2 milliseconds is a lower bound. Since the quark star 
forms from the hot neutron star after a second collapse, its ro- 
tation accelerates by a fact or ~ 1.5, according to the moments 
of inertia calculated by Bei ger & Hae nsel (2002). Thus, = 3 
ms would be a representative value of its spin period. The rota- 
tional energy of the star being W, = /*Dj/2, the kinetic energy 
available from the differential rotation is Wd = aoW,. This is 
the energy of the torsional oscillation, if it is global. For ao - 
0.025, f , = 3 ms, and 745 = 1, W/, ~ 6 10'*'' ergs. 

The oscillating toroidal field in the star is at its maximum 
amplitude Bj- when the energy of the torsional oscillation Wd 



is entirely in magnetic form, which implies that, for a global 
oscillation of a star of volume V, VB^ - ^ttWd- Thus : 



Bj 



10' 



10 a D hs 
PUms)RlJ 



gauss ■ 



(8) 



For Pt - 3 ms, ao - 0.025, and the other parameters equal to 
unity, Bt ~ 1.7x 10'^ gauss. Thus, for a poloidal magnetic field 
of order of a few lO''* gauss, the torsional Alfvenic oscillation 
is non-linear. Equation (O indicates the typical amplitude of the 
toroidal field inside the star. The field distribution has however 
a certain profile with depth and latitude, so that Eq.® is not a 
precise estimate of the sub-surface toroidal field. 

The matter's angular velocity in the star is -t-5Q, where SO. 
varies with position and time and has a null time average value. 
Indicating time averaging by brackets, we have 



W. + Wi 



(9) 



From Wd - aoW*, it is found that the r.m.s. value of the time- 
varying angular velocity equals -\fao Q* and the velocity ampli- 
tude of the torsional wave close to the surface is about a factor 
of y/oD times the rotational velocity: 



(10) 



When magnetic flux is present only in a superficial layer of mass 
m, only that part of the kinetic energy of differential rotation that 
develops in this layer feeds the energy of the torsional oscilla- 
tion. If this energy is distributed proportionally to mass, the es- 
timate given by Eq.lfTOt of the oscillation's amplitude remains 
valid. The radial component of the current in the star is then: 



c c d 

j^- — e,.- rotB - — —I sinflB,* 

47r 47rr sin 6 



^(sin.B,) 



(11) 



If this component of the current does not vanish, a DC current 
could flow from the star to the magnetospheric lepton plasma. In 
Sect. (15.31 ). we consider the consequences of these currents. 



3. Energy-extraction mechanisms 

The energy of the torsional oscillation could leak out of the star 
by a number of different mechanisms. We consider these mecha- 
nisms in turn and discuss whether they could represent the origin 
of long GRBs. 



3.1. Viscous damping of tlie torsionai osciiiation 

The oscillation could be damped in the star by viscous friction 
or Ohmic dissipation and then escape by means of heat conduc- 
tion and thermal radiation from the surface. The surface of a 
quark star at temperatures T ~ 10'" °K is an efficient s ource of y 
photons and e^e" pairs (see e.g.lAksenov et al. l (l2003 l) and refer- 
ences therein). However, Ohmic dissipation is negligible under 
the conditions assumed in Sect. ( 12.1b . The shear viscosity of nor- 
mal q uark matter has been calculated by iHeiselberg & Pethickl 
(Il993l) . The viscosity of quark matter with unpaired components 
is of a comparable order of magnitude. The Reynolds number 
for a scale R and a velocity dv^ given by Eq. ( fTOl i is found to be 

of order lO''*. Thus, shear-viscous dissipation is negligible^ 

Even superfluid quark-matter suffers bulk viscosity ( MadsenI 

II992h . An Alfven wave, being non-compressive, is not damped 
by bulk viscosity at the linear approximation, but the Alfvenic 
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torsional oscillation is non-linear. By non-resonant coupling, its 
magnetic -pressure gradients generate a compressive oscillation. 
Bulk viscosity acting on this compressive part of the non-linear 
oscillation causes a damping which can be calculated by solving 
the MHD equations perturbatively to second order We simpli- 
fied the calculation of this damping by considering an homoge- 
neous medium of mass density po contained in a Cartesian box 
with unperturbed density po and magnetic field Bq e, and an ex- 
tension L - njk in the z-direction. The solution to first order is 
the non-compressive standing Alfven wave. We define caq and 
csQ » cao to be the Alfven and the sound speed in the unper- 
turbed medium respectively, Bj the toroidal magnetic amplitude 
of the wave and ^ the coefficient of bulk viscosity. The second- 
order solution brings in the following damping time: 



Tbv - 



32po 



4 

Bj 



(12) 



The bulk viscosity in quark-matter depends on the finite time 
required by quarks to return to the weak-interaction equilib- 
rium after the flavour equilibrium is disturbed by the leptonless 
strangeness-changing reaction 



u + d s + u 



Any compression of the medium causes such a disturbance, be- 
cause the s quark is more massive than the u and d quarks. For 
colour-superconducting quark matter, the existence of a gap A 
drastically reduces the reaction rate when UbT <k A. The gap en- 
ergy A, which is uncertain, is between 1 and 50 MeV (Madsgn| 
12000 ). According to iMadsenl (il992l) . the bulk viscosity coeffi- 
cient ^ of quark matter experiencing an harmonic density pertur- 
bation of pulsation o) is: 



(14) 



where T is the temperature of the medium. Equation ( fT4] i is valid 
only when the Fermi energies of the s and d quarks differ by less 
than InksT. The coefficients a and /3 are: 

a = {4n^/9)KFmykliul , (15) 

= (7rV9) Kj klfi'', h'c' (l + my/(4ixl)f ■ (16) 

In Eqs. JTSll-dTSIl, h is the Planck's constant, c the speed of light, 
ks the Boltzmann's constant, m, the mass of the strange quark, 
and fid the Fermi energy of the d quarks, which in a IM© quark 
star of radius 10 km is 196 MeV. The mass nis is expressed in 
units of 100 MeV by 



nisC - m.5ioo X 100 MeV 



(17) 



The rate of the reactions represented by Eq. ( fTsT l depends on the 
weak-coupling parameter (Mad sen.l992i) : 



7^^ = 3.0781 lO^'^g-^cm-i^i^ 



The time needed to re-establish the equality of the chemical po- 
tentials of the s and d quarks after a perturbation (the strangeness 
equilibration time) is r,,, = {pT^Y^I^. Writing T = 10''°K, the 
strangeness equilibration time is: 



: 6.9 IQ-^T-^ i 



(19) 



This is far shorter than the period of the torsional wave, which 
means that s and d quarks always remain close to the equilib- 
rium of the reaction given by Eq. ( fT3] l when the quarks are in a 



normal state. For these representative numbers, (//, - jiafllnkBT 
remains small. This legitimates Eq. ( fT4l i. in which uP' can also 
be neglected so that: 



( = 2.45 10^^^ m.fio„ T, 



-2 -1-1 

gem s 



(20) 



If the quarks are in a colour-superc onducting state with a gap 
A, Kf is reduced by expi-2A/kBT) (lMadsenll2000h . For A = 1 
MeV and T > 2.2 10^°K, the relaxation time of reactions ( fT3] ) 
remains less than the period of the torsional wave. For lower 
temperatures, it rapidly becomes much longer For example, at 
10^ °K with a gap of 1 MeV, this time is 8 10^ sec, which is so 
long compared to the period of the wave that bulk viscosity is 
quenched. The damping time r^,,. in a 10'"* gauss poloidal field 
with a toroidal magnetic field amplitude of 10"' gauss is given 
by Eq. ( fT2b . For normal quarks this time is: 



n, = 4.8 10'" m;fo„ s 



(21) 



This is much longer than the torsional wave period, owing to 
the fact that the compression in this wave is small, so that bulk- 
viscous damping is negligible. 



(13) 3.2. Flux emergence by magnetic buoyancy 



The internal magnetic field can emerge through the surface of 
the star and expand into the quasi-vacuum outside as an elec- 
tromagnetic signal. This point o f view is adopted in the model s 
of a number of authors such as iKluzniak & Ruderman (19981), 
iRuderman etall (l2000l) . ISpruitI il999l) . and Dai & Lu ( 1998^: 
amplified magnetic fields, supposedly of order 10'^ gauss, would 
be brought to the surface of a neutron star by buoyancy forces 
and the energy will be emitted in the form of bursts and Poynting 
flux. 

We estimate the flux emergence time for a given magnetic 
field, assuming that nothing opposes buoyancy. Since the field 
in the wound-up flux tubes is essentially toroidal, these tubes 
may be regarded as thin circular annuli centred on the axis. The 
rapid rotation of the star inhibits their expansion or contraction 
perpendicular to the axis, so that the flux tubes move parallel 
to it towards the closest pole. Their length / remains constant. 
Consider a flux tube of a small cross section S and length /, 
threaded by a field B. Under perfect MHD conditions, it con- 
serves its magnetic flux and its baryonic content. For subsonic 
motions, it also remains in pressure equilibrium with its envi- 
ronment. We define Pextiz) to be the material pressure, Bextiz) to 
be the magnetic field in this environment at an altitude z along 
the polar axis, and Pi„{z) to be the material pressure and B,„(z) 
the magnetic field in the tube when it reaches the altitude z. Total 
pressure equilibrium implies that: 



PexAz) + 



87r 



Pin(z) + 



Bjiz) 
87r 



(22) 



(18) The pressure Pc in the inner regions of a quark star is about 10 



,35 



erg cm" . Since the magnetic pressure of a field of 10 gauss 
is far lower, the difference between the matter densities p,„ and 
Pext in the tube and in its environment can be calculated pertur- 
batively: 



pextiz) - Piniz) 



3(BI(z)-BUz)) 
8;rc2 



(23) 



where we have used dP/dp = c^/3. The difference in magnetic- 
energy density should be added to derive the difference in total 
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energy density. Denoting by (-gz) the component of gravity par- 
allel to the rotation axis, the vertical motion of the flux tube is 
described by the equation 



\pm(z) 



BUz) - Bjiz) 
4nc^ 



(24) 



The existence of the buoyancy instability depends on the dis- 
tribution of the magnetic field in the star. If the field intensity 
increases with altitude such that Bi„{z) is always less than Bgxtiz) 
at a higher level, the distribution of flux is stable. By contrast, if 
the flux tube moves in an unmagnetized environment, buoyancy 
can only be inhibited by density differences between the mag- 
netized and non-magnetized medium. These diff'erences may re- 
sult from a situation of non-equilibrium of weak interactions in 
the moving fluid, as described below. For Bgi/ = 0, if we as- 
sume that Bin remains almost constant during the motion and 
that ~ 0.5 (GM/R^) and p,„ ~ p. = M/V, we find from Eq. 
( l24l i that the buoyancy time Tt, needed to raise the tube by about 
a stellar radius is: 



GB^ 



2.8 10"^ s 



10'^ G 

Bin 



(25) 



This is the time required to bring an isolated flux tube to the 
surface when the toroidal field has reached the value B,„. After 
the sudden formation of the quark star, the differential rota- 
tion causes this toroidal field to increase in strength as B{t) - 
Bj sin(f/f r), where Bj- is the maximum toroidal field developed 
in the torsional oscillation represented by Eq. (O. Buoyancy 
starts to be eff'ective only when the growing toroidal field has 
reached a value such that its buoyancy time in Eq. ( 1251 ) calcu- 
lated for B(t), has become shorter than the age t of the new-born 
quark star For standard stellar parameters and Bpi4 = 1, Pt - 
3 ms and ao = 0.025, this occurs at the buoyancy starting time, 
fj, ~ 0.5 s. 

The ratio of the buoyancy time ti, (Eq. (|25T l) to the period of 
the torsional oscillation Pt (Eq. (|7J) is: 



Pt 



1 



Bp_ /GM\-i/2 

7^1' \'r?) 



(26) 



where Bp is the poloidal field and B,„ is the total field. When 
the wound-up magnetic field strength is close to its maximum 
value (Eq.®), B,„ x 100 Bp. This implies that when nothing 
opposes buoyancy, wound-up fields in isolated flux tubes float to 
the star's surface in a time t/, shorter than the period Pj of the 
oscillation. 

However, buoyancy motions are reduced or quenched when 
the ascending magnetized matter becomes denser, at the same 
total pressure, than matter in its neighbourhood. This may hap- 
pen when the magnetic field pervades the entire volume of the 
star and the field intensity increases with the altitude z- Another 
effect opposing buoyancy is when reactions such as Eq. ( fT3] ) or 
the j0-decay reactions 



d <r^ U + I 



(27) 



cannot reach equilibrium in the buoyancy time Tt,. The relax- 
ation time Tsg of the strangeness-changing reaction ( fT3] l in nor- 
mal quark matter is about 7x 10"^ Tj^ seconds (Sect. [TT). For 
colour-superconducting quark-matter with a gap energy A, the 
time to achieve equilibrium o f the same rea ctions is length- 
ened by a factor exp(2A/kBT) ( lMadsenll2O()0b . It may become 
longer than both r^, and Pp if the gap A is sufficiently large. 



Similarly, the relaxation time of the quark yS-decay reactions 
in Eq. ( |27l ) is, for normal quarks, Tg ^ 2.7 x 10"^ T^^ seconds 
(Iwamoto 1983). For colour-superconducting quark-matter, this 
time is lengthened by a factor exp{A/kBT). 

The temperature is the parameter controlling whether chem- 
ical equilibrium of the reactions in Eqs. ( [Til l and ( |27] i can be 
achieved on a given timescale. When a proto-neutron star col- 
lapses into a quark star, an energy of about 10^^^ ergs is re- 
leased, which is reflected in the initial temperature of the new- 
born o bject of ~ 3 X 10" °K. The star is then opaque to neu- 
trinos dSteiner et al.ll200il) . It cools by emitting thermal and 
Vg's from a neutrino-sphere, thermal photons of frequency higher 
than th e plasma frequency and lepton pairs. According to lUsovl 
(1200 Ih . thermal, photon emission dominates over lepton emis- 
sion at r > 5 X 10'°°K. At 10" °K, the photon emissivity is 
barely smaller than that of the black-body. Adding neutrino and 
antineutrino thermal emission, the net effective emissivity at this 
temperature is a little less than a factor of two higher than the 
black-body emissivity. The star then cools to about 3 x 10'"°K 
in 0.5 seconds. 

Does the non-equilibrium of the strangeness-changing reac- 
tions or the y6-reactions suppress the ascent of magnetized matter 
to the surface? If quarks are in a colour-superconducting state 
with a gap A, the matter in the buoyant tube retains its original 
strangeness during its ascent if Tb < Tsg expilA/ksT), where 
Tsg is the relaxation time of the reactions in Eq. ( flj] ) in normal, 
quark matter (Eq. (fT9]l). This condition is satisfied when: 



A(MeV) > TiQ (2 logjo Tk, + 3.6 - logio Bis) ■ 



(28) 



The time Tb is definitely longer than Tsg at 3 x 10 K, the tem- 
perature when buoyancy starts, meaning that in the absence of a 
gap, the reactions in Eq. ( fT3] ) remain in equilibrium. At this tem- 
perature and for Bi6 - 1, the inequality of Eq. ( |28] ) holds true 
when A is larger than about 14 MeV, which is plausible because 
the gap could be as high as 50 MeV (Madsen ,2000) . For a gap 
larger than 14 MeV, the reactions in Eq. ( [13] ) will remain out of 
equilibrium during the buoyancy motion. Similarly, the y6-decay 
reactions in Eq. ( [27b remain frozen during buoyancy motions if 
Tb < Tp expiA/ksT). This inequality is satisfied without the need 
for a gap when Tio < B^^^, which is not quite the case for T = 
3 10'"°K and Bi6 = 1; this implies that, at this temperature and 
in the absence of a gap, j6-decay reactions remain more or less 
in equilibrium during buoyancy. In the presence of a gap, the 
inequality t/, < T/sexpiA/ksT) is satisfied if the gap A is such 
that : 



A(MeV) > 2 Tio (41ogio Tio - logn, Bie) 



(29) 



At 3 X 10'"°K and for Bi6 = 1, the inequaUty ([29]) holds true 
when A is larger than 12 MeV. For this or a larger gap, the re- 
actions in Eq. ( [27] ) keep out of equilibrium during the buoyancy 
motion. It is surprizing that the gap values which freeze the re- 
actions in Eqs. (T3[ and ( [27] ) on the timescale Tb are so close. 
This results from the fact that the relaxation time of reactions in 
Eq. (T3[ lengthens more rapidly with gap energy than that for 
the y6-decay reactions in Eq. ( [27] l. 

We therefore have two situations. The gap is either less than 
10 MeV and both reactions in Eqs. ( [T3] ) and ( [27] i reach equilib- 
rium on a timescale shorter than the buoyancy timescale when 
buoyancy starts. In this case, chemical non-equilibrium has no 
role in limiting buoyancy. Otherwise, the gap exceeds 14 MeV 
and both reactions remain frozen on the buoyancy timescale. We 
disregard any intermediate situation. 
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Depressurized, non-equilibrated, quark matter weighs more 
than equilibrated matter at the same pressure because its energy 
density is not minimal. For a gap larger than 14 MeV, buoyancy 
is quenched when the total energy density e,„ in the rising flux 
tube (including its magnetic energy density) exceeds the total 
energy density eext in the ambient medium. To illustrate this, we 
consider a flux tube reaching a region in the star, at an altitude 
of z, where its material pressure is less than at the altitude zi 
where it started its ascent. We assume that during its motion the 
weak reactions remain frozen. The difference between the mass 
density p^g of equilibrated matter at this pressure and the mass 
de nsity pfr of the f rozen matter at the same pressure is expressed 
bv lHaensel & Zdunik (.2007,) as; 



Pfr - Peq - fah Peq 



(30) 



They calculated fat for cold quark-matter in which only the fi- 
decay reactions dZTb are frozen. However we are interested in 
a situation where the reactions in Eqs. ( fTlT l and ( |27] ) are both 
frozen. We find that in this case: 



=6.4 X 10-^71" 



.si 00 



(31) 



where mjioo is defined by Eq. ( fTTl ). The field B,„ of a flux tube 
that is still buoyant at an altitude z where the ambient magnetic 
field is Bext must be such that 



Bl>Bl„ + AnC^fabPeg{Pin), 



(32) 



where P,„ is the material pressure in the rising flux tube at this 
altitude. Since fab <k 1, we can assume P,„ to be almost equal 
to the total external pressure. For isolated flux tubes moving 
through an unmagnetized medium, Eq. ( [32l ) becomes: 



B]^ > AnC^fabPext: 



(33) 



where pext is the mass density in the equilibrated unmagnetized 
environment. At the star's surface p^xt - 4S, 3 being the bag 
constant. We define S50 = S/ (50 MeV fm""*). We assume that 
Ms =100 MeV. If quark matter is to be more stable than nucle- 
onic matter, S50 should not exceed 1 .6 for free quarks confined 
in the bag, and 1 .4 for the QCD coupling constant 0.2 (see, e.g.. 
Fig. 8.2 in lHaensel et"an ( l2007h ). With f^b given by Eq. (EB, the 
minimum value of a field that would be buoyant near the surface, 
at a level z, is: 



B, 



mm, top 



■■ 4 X 10'^ m'^,l S'l^ gauss ■ 



(34) 



Deeper inside the star, at a level z\, the field in this same 
flux tube had a value B„„„ due to the conservation of its mag- 
netic flux and quark number content. From Eq. (l22l l. Z?„„„ = 
Bmin,top (Pextizi) / Pextiz))^^'^ ■ The prcssure deep inside the star is 
taken to be p,c^/3 a; Mc^/AnR^. The pressure close to the star's 
surface is the bag constant S. The minimum buoyant field deep 
inside the star, B^m, is then: 



/i,„,„- 8.5 X 10 m„o||/<,n | — 



3/4 



gauss ■ 



(35) 



Taking the moment of inertia of the star to be /» = 0.4 MR^, the 
ratio of the toroidal field Bj generated by the torsional oscilla- 
tion (Eq. (O) to B„i„ is: 



Bj- 
Bniin 



0.95 



(gfl/0.1)'^ ^ /M( 



i^r R 



10 '"jlOO 



jl/4 
'50 



(36) 



Buoyancy is quenched when BTlB„,i„ < 1. For P, = 3 ms, an - 
0.025, M - 1 Mq, and S50 = 4/5, this occurs (provided that the 
colour-superconductivity gap exceeds 14 MeV) when 



m,c^ > 38 RxQ MeV 



(37) 



An important question is whether the buoyant flux is entirely ex- 
pelled out of the star with the leptonic wind or whether, although 
the magnetic field partly emerges, it remains rooted in the sub- 
surface layers. 

This depends on how rapidly the magnetic field can diff'use 
through quark matter, which itself depends on its electrical con- 
ductivity (Tg and the gradient lengthscale Im of the field. Since 
the magnetic field decreases in the flux tube during its ascent, 
its cross section at the surface cannot be smaller than when it 
started. Since rapid rotation prevents radial motions perpendic- 
ular to the rotation axis, the field scale length Im of buoyant 
flux tubes cannot diminish. The fact that the magnetic diffusion 
timescale is about 10" s for conducting quark matter implies 
that during the first few minutes after the formation of the strange 
star, the magnetic flux emerging through the star's surface as a 
result of buoyancy remains rooted in quark matter at starspots. 
In this case, magnetic activity from the torsional oscillations, as 
described below, continues after the first burst of magnetic buoy- 
ancy has brought the inner magnetic field closer to the surface. 

If matter is magnetized in bulk, buoyancy assumes the form 
of a convective instability. When it develops, the more magne- 
tized material is brought to the star's surface, while the less mag- 
netized material sinks deeper into the star. This results in a re- 
distribution of magnetic field in the star, not in a net loss of flux. 
The end result of the field redistribution should be close to a 
state of marginal buoyancy instability. A fraction of the surface 
magnetic-flux tubes should emerge from the star, baryonic mat- 
ter draining down along the field as it emerges. However, since 
this matter cannot diffuse out of the field, the emerging magnetic 
loops remain connected to the subsurface flux. The magnetized 
volume experiences little change in this process, so that a sub- 
stantial part of the star's volume remains magnetized, if it was 
initially, and magnetic activity from the torsional oscillation per- 
sists after the flux redistribution. 

In the following, we consider cases when the initial field 
stratification in the star is stable against buoyancy or when the 
colour-superconductivity gap is larger than 14 MeV and the 
mass of the strange quark sufficiently high to inhibit the buoy- 
ancy of wound-up magnetic fields. Our results also apply to 
when the star was initially magnetized throughout a substantial 
fraction of its volume and remained so after a short, first episode 
of buoyancy. 

3.3. Direct magnetic dipole radiation 

The time-dependent, internal, stellar magnetic field may be a 
source of electromagnetic emission from the star's environment, 
whether a vacuum or a leptonic plasma. For example, if the new- 
born star is an oblique rotator („ Usov..l992 ). it will emit electro- 
magnetic radiation due to the rotation of its magnetic dipole. 
We define x to be the angle between the magnetic and rota- 
tion axes. Bp the polar field, and Q, the star's rotation rate. 
The power emitted in va cuo by the magnetic dipole rotation is 
(iLandau & Lifshitzlll975h : 



(38) 



An orthogonal rotator with a polar field of lO'^ gauss and a ro- 
tation period of 3 ms would radiate a flux of 1.9 x 10"*^ erg s"'. 
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similar to that required to explain the high-energy photon emis- 
sion of the gamma-ray burst. For a polar field of only lO'"^ gauss 
the emitted power should decline to 1.9 x lO'*^ erg s which 
is insufficient to account for a GRB. We note that this emission 
taps the rotational energy of the compact star, which is about 
Eroi - 1*^1/2 » 2.2 10^' ergs for a rotation period of 3 ms. 
Even at this high rate, the radiation of a 10'^ gauss millisecond 
magnetar should last for about lO"* seconds. In the next section, 
we discuss whether the internal star's torsional oscillation could 
somehow act as a substitute for a rotating, magnetic dipole. 



4. Radiation by torsional oscillation 

The collapse leads to a state of differential rotation in the star, 
the angular velocity varying either with depth or with latitude or 
both. In an aligned rotator, somewhat analogously to the rotat- 
ing oblique dipole, the oscillating internal toroidal magnetic field 
may act as an antenna generating a large-scale electromagnetic 
wave in the star's environment at the period Pt of the torsional 
star's oscillation. We calculate in Sect. 14. li the power emitted in 
a vacuum environment. The radiation in a leptonic wind is con- 
sidered in Sect. 14.21 



4.1. Radiation in a vacuum 

To study the electromagnetic emission from the compact star 
driven by the torsional oscillation, we begin by calculating the 
electromagnetic field in an external vacuum. Maxwell's equa- 
tions are solved outside the star under the boundary conditions 
that Br and the tangential components of the electric field E are 
continuous at the star's surface. The matching of the conditions 
at the star's surface requires neither nor Bg to be continuous, 
since a surface current could support a sharp discontinuity in 
these components. We denote by a superscript < (or >) quanti- 
ties relevant to the inside (or the outside) of the star The electric 
field just below the star's surface is given by the law of perfect 
conductivity: 







(39) 



Since the velocity of the fluid in the star is assumed to be az- 
imuthal only, the condition that the tangential components of the 
electric field are continuous reduces Xo -Q and: 



cE> = -y^Bif ■ 



(40) 



In the presence of a torsional oscillation in the rotating star, 
Eq. ( |40l ) has both a time-varying and a constant component. The 
latter determines the time-independent, external electric field, 
while the former determines the outside radiation caused by the 
torsional oscillation. The boundary condition in Eq. ( l40b deter- 
mines completely the solution in the vacuum outside the star. We 
calculate the electromagnetic field radiated out of the star by the 
internal torsional oscillation of pulsation = Itz/Pt, assuming 
axisymmetry. The toroidal field is then the only time-dependent 
component of the outer magnetic field. Omitting for simplicity 
the superscripts > which refer to the outside region, the equa- 
tions for the electromagnetic fields in this region are Ampere's 
and Faraday's equations in a vacuum, given by: 



cV A B 

dt 



0. 



dB 

— H- cV A £■ = 

dt 



(41) 



(42) 



The only non-vanishing and time-dependent component of the 
magnetic field is the azimuthal one. It is useful to introduce an 
angular potential ^(r, 9, t) such that 



(43) 



The system (l4TI)-(l42b reduces to an equation for Z?^ alone, which, 
for harmonic time-dependance in the form of exp(-;Qj-f), trans- 
lates into the Helmholtz equation for //: 



(44) 



The operator A is the ordinary scalar Laplacian. We expand 
li(r, 0, t) in spherical axisymmetric harmonics. The solution for 
each harmonic component of degree { of ji is, to an arbitrary 
multiplicative factor: 



<i/2(^"'-) 



exp(-iD.T t) P"(cos 0) 



(45) 



where ko is given by the dispersion law of free-space electromag- 
netic waves, that is kg - 0.t/c. In Eq. ( l45T l. H^^^^^j(x) is a semi- 
integer Hankel function and f ° is the Legendre polynomial of or- 
der (. The value of f = 2 is the lowest value of { pertinent to our 
problem. In fact Eq. (|6]l illustrates that Z?^ inside the star is gen- 
erated by ABp terms. The dipole components of the poloidal 
magnetic field Bp correspond to the lowest value of {, {I - 1), 
since - B^ cos and Bg - -B^ sin 0, where Bp is the polar 
field. Similarly, the lowest value of f for the time-dependent rota- 
tion velocity at the star's surface is( - 1, when dV^ - Svf,^ sin 0, 
6vgg being the velocity of the time-dependent part of the rota- 
tion at the equator This angular variation in SV^ corresponds 
to a constant-amplitude modulation of the rigid-body rotation- 
rate at the star's surface: 6Q.(R, 0) - 6Q.eq. Such an oscillation 
would be induced by variations in the fluid's angular velocity Q 
with depth. Differential rotation in latitude corresponds to ^ > 2. 
Since the vector product of dV^ and the poloidal magnetic field 
is expanded in spherical harmonics with { > 2, we now restrict 
our attention to emission in the ( - 2 mode. The semi-integer 
Hankel functions can be expressed as a sum of a finite number 
of simple terms. For 1-2, the solution for which behaves as 
an outgoing wave at infinity is: 



k,i . kQ 
+ 1^ - 



(1 - 3cos2 0)e 



i(k„r-nTt) 



(46) 



where mo is a complex factor. The complete solution can be de- 
rived from Eqs. ( |4TI ). ( l42t . and ( 1431 ) and is written in the follow- 
ing form, where Bo is a complex amplitude: 



kn . 

+ - 



^2Q eHhr-ClTt) 



(47) 



= ^ ( ^1 + 4 - (1 - 3 cos^ 0) e*'-"^", (48) 



E>^ 



R^Bqc 
2iQ.T 



+ ikn 



3r 



sin 20 e 



i{k„r-aTt) 



(49) 



The relations (|47l)-(|49l) solve the system of Eqs. (I4n)-(l42t. The 
complex amplitude Bq is determined from the boundary condi- 
tions, such that the component of the electric field is contin- 
uous at the star's surface (Eq. (|40l)). The time-dependent field 



J. Heyvaerts et al.: Quark star undergoing torsional oscillations 



9 



component Eg of Eq. ( |49l ) must match the corresponding £ - 2 
time-dependent part of the field component Eg just below the 
star's surface, which requires that: 

c Eg = Bp Sv^eq sin ecos 6 , (50) 

where Bp is the polar field and dv^eq is the time-dependent part 
of the { - I component of the azimutal velocity at the equator 
By matching the two members of Eq. ( fSOl l we derive the complex 
wave amplitude By. Since the star's radius is far smaller than the 
wavelength of the emitted wave, Eq. ( |50l ) should be evaluated 
to the dominant order in the small parameter koR ^ 2 x 10"^, 
providing: 

B^JMbp'-^. (51) 
2 c 

The velocity amplitude of the torsional oscillation is given by 
Eq. ( [Tol l and k„ - Q.t/c. The modulus of Bq is then: 

I Bo I = 6.94 X 10^ ^^Jl^jj^jB^,4i;2„ gauss .(52) 

The magnetic amplitude of the wave is far smaller than the 
sub-surface magnetic field because of the significant impedance 
mismatch between the star's interior and the outside vacuum. 
Denoting by the Alfven speed inside the star, the ratio of 
these impedances is v^/c a: 4x 10^^ Bpu- The peculiarities of 
the spherical wave solution for Eg are also responsible for the 
smallness of this amplitude, which is not set by assuming conti- 
nuity of B^. The correct boundary condition is Eq. ( |50l l and its 
fulfilment implies that B^ is discontinuous at the star's surface. 

It is interesting to evaluate the Poynting power radiated off 
the star's surface. If the low-frequency wave emission can be 
represented by radiation in vacuo, the solution of Eqs. (|47Ti-(|49]l 
provides an upper bound to the power that may be dissipated in 
the star's environment and radiated away as X and y photons. 
The radial component OJ, of the Poynting vector associated with 
the low-frequency radiation is 

<i>'p^^[^'Re{E;B^)y (53) 

where the superscript * designates the complex conjugate and Ke 
the real part of a complex number. The power Tyac radiated as 
Poynting flux by the torsional oscillation in its supposedly vac- 
uum environment can be calculated from Eq. ( l53b by integrating 
cDp over the star's surface, taking Eqs. (I47b-(l49b into account. A 
number of simplifications occur in this calculation, which finally 
infers that: 

'P,ac^-^[\Bo\R'^lf. (54) 

The magnetic amplitude Bq of the wave is given by Eq. ( fSTT l and 
6v0eg is given by Eq. ( fTOb . We then have: 

^™-=5i^^^"'"*^" ^''^ 
Numerically, the power !Pi,3c amounts to: 

= 2.06 X 10- Bl,,R- (^) (^f (^f erg s^ • (56) 

A glance at Eqs. ( [38] l and (|55l l indicates that much less energy 
is radiated away in an outside vacuum by the torsional oscilla- 
tion than by an oblique, rotating, magnetic dipole. The Poynting 



power radiated by the torsional wave is smaller than the emis- 
sion of a rotating dipole for several reasons. First, the radiation 
is quadrupolar instead of dipolar. Then, the field | Bq I is only 
of the order 10*^ gauss for the values of parameters adopted as 
representative (Eq. (l52b). much less than the polar field of an or- 
dinary pulsar or magnetar Finally, the period Pj- is of the order 
of a few seconds, much longer than the rotation period of the 
new-born compact star As a result, the power in Eq. ( |56] | falls 
short by many orders of magnitude of the observed power of y 
radiation in a long GRB ! 

4.2. Radiation in a leptonic wind 

Could the presence of a circumstellar, leptonic plasma drasti- 
cally change the power radiated by the torsional oscillation? 
This plasma originates in charges, electrons, and positrons that 
have passed the bag of the quark star We note that at a dis- 
tance from the star of larger than the light-cylinder radius, 
the pl asma cannot be in rigid corotation, but must flow out- 
ward (IGoldreich & Julianl ll969). Since the wavelength associ- 
ated with the frequency of the torsional oscillation is much larger 
than the light-cylinder radius of the rapidly spinning star, the 
wave emitted by this oscillation propagates into the wind driven 
by the rapid global rotation. We have to determine the amount 
of energy of the torsional oscillation radiated per second in these 
conditions. We assume that the wind has already reached its ter- 
minal velocity at the surface of a sphere of radius comparable 
to that of the light-cylinder, c/Q*. Since the ratio of the torsional 
oscillation period to the spin rotation period is large, the torsional 
oscillation appears, on both of the scales of the spin period and 
the light-cylinder radius, as a quasistatic perturbation. Its effect 
is not only to emit a signal that assumes the character of a wave 
at distances larger than c/Qt, but it also modulates the wind in 
which it propagates as a result of the variations imposed on the 
conditions of its lauching. These modulational effects are dis- 
tinct from emission of low-frequency radiation in a given wind. 
Much of the action causing wind modulation occurs below or 
close to the light-cylinder and will be discussed in Sect. |5] 

We now calculate the emission by the torsional oscillation in 
an expanding, possibly resistive, leptonic wind. The conductiv- 
ity cr of the medium is assumed to be real. This is because the 
wave would be highly non-linear if the gyrofrequency of lep- 
tons in the wave's magnetic field was much larger than the wave 
frequency. As a result, the effect of the plasma current on the 
real part of the index of refraction would become negligible and 
the wave wou ld force its way non-linearly th rough the leptonic 
environment ( Asseo et al.lll975l:ISalvatiiri978h . We may then re- 
strict our consideration to resistive effects. The modulus of the 
wind speed w is assumed to be constant, both in time and space, 
and oriented radially outwards: w = wCr. We assume the wind 
to be ultra-relativistic and to have a velocity equal to the speed 
of light. When taking the wind to be radial, we assume that coro- 
tation is lost at distances of the order or larger than c/O^. The 
background, magnetic field in the wind is severely wound up 
by the star's rotation. At distances much larger than the light- 
cylinder radius, its azimuthal component dominates over the 
poloidal component and declines proportionally to l/r. We thus 
neglect the poloidal field component and assume the unperturbed 
magnetic field to be azimuthal, so that: 

Bo^Bo-me^ ■ (57) 
r 

Specifically, we consider /(0) - cos0/{\ cos 6 | sinff). The elec- 
tric current associated with Eq. ( |57] i is radial. For the adopted an- 
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gular profile it reduces to zero almost ever ywhere, except at both 
the p olar axis and the equatorial plane. i Hevvaerts & Normaiil 
(|2003) demonstrated that the magnetic field in a perfect MHD 
wind asymptotically becomes potential almost everywhere, the 
cuiTent being confined to boundary layers about the polar axis 
and at surfaces where the poloidal polarity reverses. Our choice 
of f{Q) complies with this. The singularity at the polar axis rep- 
resents the current carried by a jet, while the change in the di- 
rection of the field at the crossing of the equator is caused by the 
change in polarity of the poloidal field at = 7r/2. Ohm's law 
infers that: 



j - a{E +-whB^+ peW 



(58) 



By considering the divergence of Eq. (I58l l and solving the re- 
sulting differential equation for the charge density p^., it is found 
that the latter vanishes in the unperturbed wind. The unperturbed 
electric field then vanishes too. We define E, B, j, p^, and u to 
be the perturbations of the electric and magnetic field, current 
density, charge density, and leptonic fluid velocity, respectively. 
It is sufficient to describe the lepton's dynamics at the inertia- 
less (also called force-free) approximation. The magnetic-field 
perturbation is toroidal in the considered geometry and the elec- 
tric force peE is a negligible second-order term. The perturbed 
Maxwell equations. Ohm's law, and the dynamical equation in 
the zero-inertia limit can be written (using the compact notation 
dt for time derivatives) as: 



cV A £■ = -d,B , 
cV A B = 47Tj + d,E , 

V ■ £■ = 47TP, , 

cj — a-{cE + w A B + u A B()j + cp^w . 
jABo = ■ 



(59) 
(60) 
(61) 
(62) 

(63) 
(64) 



The toroidal components of Eqs. f5% and ( |63] ) imply that the 
electric field is only poloidal. The other two components of Eq. 
( |63] ) infer the fluid velocity once the solution for the other un- 
knowns has been found. The current j is eliminated by taking 
the vector product of Eq. ( |60l l and Bq. As a result the conductiv- 
ity is eliminated from the equations describing the perturbation. 
With Eq. ( |59] l, this infers an equation for the magnetic perturba- 
tion, which is expressed most accurately in terms of the angular 
potential fi of this perturbation (Eq. ( |43]) ). For harmonic time- 
dependence exp(-/f2rf), it is found that ju satisfies the Helmholtz 
equation (Eq. (l44b). Regardless of the conductivity, the perturba- 
tions propagate in this geometry as electromagnetic waves, pro- 
vided the inertia-less limit is considered. The complete solution 
is identical to the vacuum result (Eqs. (|47]i-(|49]l) as is of course 
the Poynting flux at the star's surface and the emitted Poynting 
power (Eq. (l55Tl). 

5. Modulation by the torsional oscillation of the 
energy emitted in the rotator's wind 

A rapidly-spinning aligned rotator emits a wind carrying power 
in electromagnetic, potential, thermal and kinetic energy form. 
The contributions of these different forms of energy depend on 
the distance to the star Some forms of energy may dissipate en 



route or at terminal shocks, producing observable X and y radi- 
ation. Close to the compact star, much of this flux is in Poynting 
form because the kinetic energy remains low while the wind has 
not yet been effectively accelerated. The thermal and gravita- 
tional energy fluxes often constitute but a little part of the total 
energy flux. The energy output of the star in its wind environ- 
ment then enters the latter as DC Poynting flux, the radial com- 
ponent of which is given in terms of the field components just 
above the star's surface by: 

<lynr^—E>B;- (65) 



4n 



The boundary condition at the star's surface implies that Eg = 
Eg. We define to be the subsurface fluid velocity and B, the 
radial field component, which is continuous across the star's sur- 
face. We then find that: 

Ooc = - ^ MR) BAR) b;(R) ■ (66) 

The value of B^ at the base of a relativistic wind is given approxi- 
mately by Eq. ( |67] l below. This can be sketchily explained as fol- 
lows: due to the rotation of the star and the effect of flux freezing, 
a toroidal field is generated on open field lines from the poloidal 
field. The knowledge that the foot point of a field line is rotating 
is propagated along this line at a finite velocity Vprop, by means of 
convective transport and propagation as an Alfvenic signal. The 
field line thus curves away from the sense of rotation at an an- 
gle of ^ to the radial direction, such that tan ^ = 0,7? sin 0/ Vprop ■ 
Since the wind is relativistic and the Alfven speed in the tenuous 
external plasma is close to the speed of light, Vprop ~ c and the 
toroidal field just above the star's surface is: 



R Q, sin e 
b;(R) * BAR) 



(67) 



A more precise justification of the approximate relation in Eq. 
(l67T i is omitted for conciseness. Then, from Eq. ( l66] i: 

Oflc * -T^ ■ (68) 

4-7TC 

This flux is only emitted from the polar caps, the regions on the 
stellar surface connected to open, field lines. The magnetosphere 
is closed, where the apex of the local field line is at a distance D 
smaller than the light-cylinder radius c/Q,. When field lines are 
dipolar, the polar caps extend to a colatitude 9pc, which we refer 
to as classical, and described by: 

Q.,R 



sin^e, 



pc 



(69) 



Under certain conditions however, the magnetospheric field may 
depart considerably from dipolarity (Sect. |5.3l l. When the flux 
is distributed on the star as a dipolar field, the radial, field com- 
ponent varies with 6 as Br = Bp cos 0, Bp being the field at the 
pole. By considering v^(6) to be the solid body rotation velocity 
at the angular speed Q*, we obtain, by integrating over the co- 
latitudes corresponding to the two polar caps, the DC Poynting 
power emitted by the star under these conditions: 



B^R'^Sll 
4c3 



(70) 



The power represented by Eq. dTOl l is comparable to the power 
emitted by an oblique, rotat ing dipole (Eg . (|38] |). This is a clas- 
sical result (see for example iMichell (1199 ll) ). With a rotation pe- 
riod of 3 ms, the power Vsoi ~ 1.78 x W^^B^^^^ erg s"'. This is 
insufficient to match the high lumino sity of a GR B unless fields 
in excess of 10 gauss are involved ( IUsovlll992l) . Could diff^er- 
ential rotation drastically change this result ? 
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5.1. Quasistatic moduiation of the wind 

We assume that the star experiences a torsional oscillation. The 
velocity v^{0) then differs from the solid-body rotation veloc- 
ity and varies with time. Because the period of the oscillation is 
much longer than the mean spin period, this causes a quasistatic 
change in both the structure of the magnetosphere and the po- 
lar cap angle. It even causes a change, which we neglect, in the 
shape of the light cylinder The structure of the magnetosphere 
and the energy output of the wind adjust to equilibrium values 
corresponding to the instantaneous velocity profile on the star's 
surface. 

This profile may be even with respect to the equator, or 
odd, or a mixture of both. An even oscillation is one in which 
the azimuthal velocity perturbation is symmetric with respect to 
the equator, i.e. where the time-varying azimuthal velocity is in 
phase at two points positioned symmetrically with respect to the 
equator. An odd oscillation is one in which it is antisymmetri- 
cal, i.e. where the velocity is in phase-opposition at two points 
symmetric with respect to the equator. It may appear that na- 
ture should provide only even profiles, by a principle, or rather a 
postulate, of north-south symmetry. However, this symmetry is 
broken in the case of the collapse of supernovae. It is indeed well 
known that new-born neutron stars receive a kick, that is, a net 
thrust from the collapse. This would be forbidden by the prin- 
ciple of north-south symmetry. Therefore, it cannot be excluded 
that, similarly, odd modes of torsional oscillation are present in 
the initial excitation of a new-born compact star. The amplitude 
of these modes is expected to be small, but we show below that 
it need not be large to produce important effects. The kick re- 
ceived by a new-born neutron star produces a velocity of the 
order of 200 - 500 km s"' . This corresponds to an asymmetry in 
the momentum emission of the order of IMqX 300 km s"' . The 
supernova explosion emits a momentum per steradian of about 
10 MqX 30 000 km s"'/ 4 tt. The asymmetry in the momentum 
emission appears to be a fraction of between a few 10""^ and a 
few lO""* of the total. A similar fraction of the total star's rota- 
tional energy may appear after the collapse in the form of odd 
differential rotation. 



5.2. Moduiation of the wind by an even osciiiation 

An even torsional oscillation has but little effect on the structure 
of the magnetosphere and wind because the two footpoints of a 
closed field line follow the same motion exactly if, as assumed 
in this subsection, the poloidal field lines are strictly symmetric 
with respect to the equator Otherwise, these field lines would 
undergo a twist in the presence of an even torsional oscillation, 
because their footpoints would not be at exactly opposite lati- 
tudes, and would be carried in the azimuthal direction at different 
angular velocities. It is difficult to realistically anticipate the de- 
gree af asymmetry in poloidal field lines. It could vary from very 
little to a complete absence of symmetry. The degree of asym- 
metry necessary for the magnetosphere to open will be estimated 
in Sect. (|5.3l l. Strictly symmetric field lines whose footpoints are 
moved by an even torsional oscillation are, however, not twisted 
and, as a result, no poloidal electric current is driven in the mag- 
netosphere. Nevertheless, because the rotation rate on the star 
varies with colatitude, the Goldreich-Julian charge distribution 
in the magnetosphere differs slightly from the case of solid body 
rotation, as well as the DC Poynting flux (Eq.(|68]l). Using as a 
model the following even differential rotation: 



we calculate from Eq. ( |68T l and integrate over the classical 
polar caps to derive the emitted power Peven- 

"Peven = sol (l + COS Oj-f + S2 COS^ i^r?) , (72) 

where 

5 V c \ Q. / ' " 3 \ c /\ O." 



^.^J^f^l s.-U^M'-S^] (73) 



The wind power is modulated slightly at a level of ei, which is 
about 10 % for a period of 3 ms and <5Q,/Qt = -sfaD ~ 1/3. 
The change in the time-averaged power is at the 0.1% level for 
the same figures. These small changes cannot account for the 
existence of a gamma-ray burst. 

5.3. Magnetosphere opening by an odd osciiiation 

An odd oscillation differs from an even one in that the two foot- 
points of a closed field line experience differential motion in 
longitude, introducing a twist in this field line. This causes a 
poloidal current to flow in the closed magnetosphere and drasti- 
cally changes its structure. When the twist exceeds a threshold of 
order n, the magnetosphere opens. Field expansion by the shear- 
ing of the footpoints of field lines w as first discusse d in the con - 
text of solar flares (HeyvaertsitaD[i982; Aly 1985; Lowlll990l) . 
It was established that it occurs in Cartesian geometry with a 
direction of invariance by th eorems constraini ng the properties 
of line-tied for ce-free fields (lAlvlll985L Il990l) and by numeri- 
cal simulations dBiskamp & WelterllI98 9*). The same process has 
been considered also in the case of axisymmetric structures ex- 
tending above a spherical surface on which the field lines are 
tied. In general, it was demons trated that rapid expansion occurs 
when a finite shear is reac hed (|Alvlll995j). This is also supported 
by numerical simulations (Mikic & Linkeij[T994 ). Full opening 
occurs for a finite tw ist (of order n) in som e specific examples 
dLvnden-BeU & Boilvil l994lWolfsonlll995l) . In the present con- 
text, the light-cylinder radius imposes a limit on the distance to 
the apex of closed field lines, such that field opening is even 
easier when magnetospheric inflation proceeds. When the field 
opening becomes almost complete, it causes a growth in the po- 
lar caps and the emitted wind power. In this case, the role of 
the torsional oscillation is not to modulate the energy output by 
the addition of its own electromagnetic emission but to open the 
door for a more significant wind emission from the central ob- 
ject. This enhanced wind emission would acquire its energy di- 
rectly from the rotational kinetic energy of the star, not only from 
the energy of the differential rotation, and lasts for as long as the 
torsional oscillation survives with sufficient amplitude. 

The magnetosphere opens if the difference in longitude t// be- 
tween the two conjugate footpoints of a field line exceeds typi- 
cally half a turn (i// > n). To justify this statement, one should try 
to solve for the structure of the magnetosphere as a function of 
the difference in longitude between the footpoints of field lines. 
By assuming axisymmetry, the poloidal field is represented by a 
flux function a(r, 0), so that the total magnetic field can be writ- 
ten as: 



Vfl A Cj, 

r sm 



(74) 



7? sin e (O, + sin 6 cos Qj-f) 



(71) 



Any field line follows a surface of constant a (a magnetic sur- 
face) because the magnetic flux being transmitted through a 
circle perpendicular to and centred on the polar axis passing 
at (r, 0), is (I'm = 2na{r,ff). The perturbations of the closed 
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magnetosphere of the star are of quasi-static nature (Sect. 15. It . 
Neglecting the particle's inertia, the force equation for the in- 
stantaneous equilibrium is the force-free equation; 



cpeE +j A B - 



(75) 



The electromagnetic state of the magnetosphere is not described 
by the magnetic field alone but also by the electric potential 
U(r, 6). Equation ( |75] | is supplemented by the time-independent 
Maxwell's equations; 



V A B - Anjic , 
V ■ £■ = Anpe , 
£ = - Vf/. 



(76) 
(77) 
(78) 



The components of Eq. ( |75l l can be expressed in terms of the 
functions a(r, 6), U(r, 9), and /(r, 0), the latter being defined by; 



/(r,0) = cr sin 6*50 



(79) 



The current through the circle perpendicular to and centred on 
the polar axis passing at (r, 0) is 7 = 7/2. We refer to I as the 
poloidal current. In an axisymmetric state, the electric field is 
poloidal. The toroidal component of Eq. dTSl ) then shows that 
the gradients of / and a are everywhere parallel, which implies 
that / is a function of a: I{r, 0) - I{a{r, ff)). The poloidal part of 
Eq. ( ITSI i can then be written as; 



sin^ AU VU 



where A is the scalar Laplacian and 



A*a 



sm0 d j 1 da 
50\sin0 00 



(80) 



(81) 



Eq. (IHOl l indicates that the gradients of U and a are everywhere 
parallel, which implies that f/ is a function of a, U{r, 6) - 
U{a(r, 6)). The rotation rate of the matter, Q.{d), is given by the 
electric drift velocity of particles and is found to be; 



n(a) = c U'(a) 



(82) 



The proje ction of Eq. ()80l l onto Va provides the so-called pulsar 
equation (lMichellll99lb ; 



I P-D? sin^ 6*^ 
1 ^ Afl 



2 da 
r dr 



sin^ 



2 cos da 
sin 89 



U'ia) 



(83) 



This equation has a singularity at the light-cylinder, 
which causes any field Une reaching this limit to diverge 
dContopoulos et all 1 19991) . To determine approximately which 
field lines become open, it suffices to solve Eq. ( l83l l for 
rQsinfl « c and confirm which field lines reach a distance 
larger than c/Q,. In this limit, Eq. ( [83] l reduces to; 



d^a sin 6* d I 

~P~'d0\lm0 89 



1 da\ U'ia) 



= ■ 



(84) 



It is shown in Appendix lAl that, for a self-similar model of the 
magnetospheric field, there is no solution to Eq.(l84ll with closed 
field lines when the twist exceeds n. 

We represent odd differential rotation by the following sim- 
ple model for surface differential rotation; 



6Q.{0) = sin26'cosQ7-f 



(85) 



The magnetic flux distribution on the surface of the star is a func- 
tion of the colatitude 9, so that is a known function of a, 0(a). 
The twist at time t associated with Eq. dSSl l is; 



ip{a) - — — sin20(fl) sinQj-f ■ 
£27- 



(86) 



It is implied here that the footpoint P2 is at the colatitude 9 and 
that Pi and Po are at the same longitude at f = 0. In the model 
presented in Appendix lAl we regard < 1/' > at time t as being 
half the maximum twist implied by Eq. ( |86] l. The limit value 
of n for the twist is reached for a rather small amplitude of the 
torsional oscillation because the period of the latter is long com- 
pared to the spin period. Similarly to Eq. (fTOl i, we parametrize 
as 6Q.t - y/Oodd ^* in terms of the fraction a„dd of the 
star's rotational energy available in this odd oscillation mode. 
The magnetosphere is in an open magnetic configuration when; 



sinQj-r > 10 



10"^ sec 



10 sec 



io-4y/2 



Oodd 



(87) 



As explained in Sect. 15. II a fraction 10 of the star's rotational 
energy may be stored in odd torsional oscillation modes. At 
this level of excitation of odd modes, the magnetosphere would 
be open during a large fraction of the oscillation period. There 
would be no opening only when Oodd is very small, i.e.; 



a odd < 10 



10" 



2 / 1 A \2 

10 sec\ 



(88) 



We assume that odd oscillation modes are initially excited to a 
level higher than the limit indicated by Eq. dSHT l. 

A similar result is obtained when an even oscillation is con- 
sidered (with an amplitude given by the larger value indicated in 
Eq. ( fTOb ) but the dipolar-like magnetic-field lines are not strictly 
symmetric with respect to the equator This would happen if, for 
example, the field is a non-centred dipole. We define S0 to be the 
difference of the absolute values of the latitudes of two conju- 
gate footpoints. The twist experienced by these footpoints will 
be larger than 7:, and thus the magnetosphere will open, when 
s/odO., (69 1 (n 12)) Pj > n, which translates into the condition; 



50 



> 1.5 X 10" 



10"^ sec 



10 sec \/ 10 



ao 



1/2 



(89) 



For Eq. ( 189b to be satisfied for typical values of P,, Pj, and ao, 
it suffices that the dipole field be decentred by a fraction of a few 
10"^ of the stellar radius. Since the odd torsional oscillation has 
an amplitude larger than indicated by Eq. ( [88] l or the magnetic 
geometry is north-south asymmetric to a degree larger than in- 
dicated by Eq. (|89] |. the magnetosphere will alternate during the 
oscillation cycle between a classical state, which we refer to as 
closed, and an open state. 

For an odd torsional oscillation, the configuration is closed 
when the condition of Eq. dSTT i is not satisfied. It becomes 
an open state, where all field lines are open and carry winds, 
when the twist is sufficiently large for Eq. dSTl i to be satisfied. 
During closed episodes, the polar caps opening is limited to 
0pc - (Q.,R/c)^^^, and during open episodes 0pc = n/2. There 
is a transitory state which we neglect because it lasts much less 
than a wave period. In an open state, the power fed by the com- 
pact star into its relativistic wind is much larger than the classical 
value given by Eq. dTOt . This may be the reason why the emitted 
power is enhanced considerably in the first moments after the 
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collapse, an enhancement that should dechne as the star's rota- 
tion decelerates and last at most until the odd mode amplitude 
has decreased below the limit fixed by Eq. dSSl l. We calculate 
the lifetime of odd oscillations and their associated emission in 
Sect l5.4l The idea that a GRB would be the result of pulsar-type 
emission from a compac t star with an entirely open magneto- 
sphere was considered bv'Ruder man et all (1200 0). who however 
regard the expansion of the magnetosphere as being caused by 
magnetic buoyancy rather than by twisting, as we suggest in this 
paper 

The power emitted at time t is calculated by integrating the 
Poynting flux given in Eq. ( |68] | over the wind-emitting star sur- 
face, where is given by (see Eq.dSSll): 



= RQ.^ sin (Q, + 60.^ sin 26* cos Qjf) 



(90) 



When the magnetoshere is in a closed state, the emitted power 
is, neglecting terms of order (5fi,/Q»)^: 



'Podd/d 



B^R^^t 
4c3 



Similarly, when the magnetosphere is open, it is given by: 



2 BlR^^l^ 



roddiop - j3 ^ 

Numerically, with Bpi4 - B^,/10''* gauss and ^ = 10 km: 



Voddiop^l.^lO"^ B 



48 d2 / 10 ^ s\ 1 



Pl4 I 



ergs 



(91) 



(92) 



(93) 



The power emitted during the open episodes is larger than that 
emitted during the closed episodes by a factor of 8c^/(150j/?^). 
For f , = 3 ms, this factor is about 110. For the same rotation 
period and Bpu - 3, Voddjop reaches 1.8 x lO'*^ erg s"', which is 
enough to explain the GRB emission, allowing for a conversion 
factor from kinetic wind energy to radiation that is smaller than 
unity. 

5.4. Damping of rotation and odd oscillation 

Neither the amplitude of the odd torsional oscillation nor the 
rapid stellar rotation will last long in the presence of such large 
losses. A rapidly spinning star with an odd oscillation is char- 
acterized by two parameters, the average star-rotation rate O, 
and the amplitude of the differential rotation Due to wind 
losses, both decrease in time. To calculate their evolution, a 
model of the internal magnetic field of the star is required. 
Although this is not an accurate model for the winding-up of 
the field when the rotation depends only on the distance to the 
axis, we shall assume for simplicity that the unperturbed, mag- 
netic field is uniform in the star and parallel to the rotation axis, 
i.e. that B< - BpC^. The normal component of this field on the 
star's surface equals that of a dipolar field and the inner field Bp 
equals the outer field at the z > pole. It is convenient for us to 
use the cylindrical coordinates D, (p, and z, the parameter r rep- 
resenting the spherical distance to the star's centre, R the radius 
of the star, and 9 the colatitude. Since the fluid is almost incom- 
pressible with a uniform mass density p, we assume as in Sect. 
|2]that its velocity is azimuthal and can be written as: 



y = (Q.D + vv(D,z,r)) 



(94) 



The velocity is supposedly odd in z- When writing the even 
part of the rotation velocity as Q»D, we neglect the even tor- 
sional modes, which play no role in the magnetosphere opening 



in the case of a north-south, symmetrical, magnetic structure. 
The magnetic field in the presence of the perturbation develops 
an azimuthal component equal to: 



B = Bp + B^(D, z, t) 



(95) 



For the assumed, uniform, unperturbed field, Eqs. (|5])-(|6l) can be 
written as: 



dv^ _ Bp dB^ 
'di ~ ~ 



dB^ 
dt 



= B 



4kp dz 



(96) 
(97) 



If were, at any given time, structured according to the 
Proudman-Taylor theorem, it would be cylindrical, that is, it 
would depend only on the distance to the rotation axis. Then, 
for a completely axial field as assumed above, the term on the 
right-hand side of Eq. WT\ would vanish. The assumption of a 
uniform axial field, however, is only applied to simplify the fol- 
lowing calculations. Nothing constrains the structure of the field 
inside the star, especially when it is not dynamically significant. 
For non-axial Bp, Eqs. (I96b-(|97]) still provide a sketchy repre- 
sentation of the torsional oscillation: the operator Bp djdz must 
be assumed to represent (Bp-V) (see Eqs. (|5]l-(l6ll) and the vari- 
ations in r sin 9 along a field line should be ignored, although 
they exist. Using Eqs. (I96]l-(l97l) is equivalent here to replacing 
a curved poloidal magnetic field in a cylindrical velocity field 
with a uniform axial field in a z-dependent velocity field. An odd 
torsional oscillation in a cylindrical velocity field would corre- 
spond to north-south asymmetric poloidal field lines. The simple 
approach implied by Eqs. (|96]|-(|97]| should be sufficient for our 
purposes of estimating the damping time of the torsional oscilla- 
tion. These equations hold for a linear as well as for a non-linear 
axisymmetric perturbation of an incompressible medium. The 
velocity perturbation satisfies the Alfven propagation equa- 
tion: 



B\ d\ 



5% 
df- " 47rp dz^ 

A standing wave solution of Eq. 
v^(D, z, t) - v^iD) sin kz cos Qpf , 



(98) 



(99) 



U„2 



where k and Q.t are related by the dispersion relation Qj = k v 
and = Bp/(4np). The wavenumber k could depend on the 
distance D to the axis, each cylindrical magnetic surface then 
having its own oscillation period. To keep things simple, we as- 
sume that A: is a constant equal to n/R, Qp = ttva/R, and is an 
odd function of z. The proper choice of v^iD) in Eq. (|99T l ensures 
that the velocity field of the perturbation at the star's surface co- 
incides with Eq. (|90] |. that is: 



v^(Rsmff) - 



2RSQ., sin^ecosfl 



sin(;r cos 9) 

Equations ( |97] i and ( |99] l infer the magnetic perturbation to be 



(100) 



V0 

B^ - Bp — cosfe sinQpf 
Va 



(101) 



To obtain the angular-momentum balance equation for both the 
z > and z < hemispheres, the torques acting on each must 
be calculated. Each hemisphere experiences volume torques ex- 
erted by magnetic tension and surface torques caused by the drag 
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produced by the emission of Poynting energy at the star's sur- 
face. The angular momentum J+ of the z > hemisphere and 
the angular momentum J of the z < hemisphere can be calcu- 
lated from Eqs. (gUl, and ( fTUOl i. They ai-e: 



J+ - I\D.^ + h 5Q.t cos Qj-f , 
J - /i - h SQ.t cos Oj-f ■ 

The moments of inertia I\ and I2 are: 
4;r , 



si 
Jo 



sin'*^ cos^0 tan 



(^cosfl) 



d6- 



(102) 
(103) 



(104) 
(105) 



The numerical value of the integral on the right-hand side of 
Eq. ( IIO5I 1 is 0.187. The moment of inertia /* of the entire star 
with respect to the rotation axis is /» - 2I\. The torque Tb+ 
exerted by magnetic tension on the z > hemisphere can be 
calculated from the Lorentz force density (Eq. (|96l)). If is 
given by Eq. (llOlb . this implies that: 



B^Ii^ r"/^ 111 \ 

— — JQ.sinQrf | sin'*6icos^6itan(-cos6l)t/0-(lO6) 
va Jo \2 / 



The magnetic tension torque Tb- exerted on the z < hemi- 
sphere is Tb- - -Tb^- The Poynting torque dTp^ exerted on 
the strip of the star's surface between colatitudes 6 and 6 + d9 
in the z > hemisphere is related to the Poynting power d'P+ 
emanating from that strip by d'P+ = - dTp+, where is the an- 
gular velocity of the fluid at that colatitude and at that time (see 
Eqs. (|90l), dMll and (l99]l-(fT00li). The torque dTp^ on the strip 
between and 6 + dO in the z < hemisphere is similarly calcu- 
lated: 



dTp+^ --l^(^l,+6£l,sm26cosQTt)sm^0cos^ed0, (107) 
2c 



- 6£l,sm26cosQTt)sm^ecos^ed0 ■ (108) 

2c 



dTp 



The total Poynting torques Tp+ and Tp- on the z > and z < 
hemispheres are derived by integrating Eqs. ( 1107b and ( 1108b . re- 
spectively, over colatitudes from zero to the polar cap angle 6pe. 
When the magnetosphere is closed. Ope x (Q,/;/c)'^^. When it is 
entirely open. Ope - jt/2. The angular-momentum balance equa- 
tion for the z > and z < hemisphere are given respectively 
by: 



dJ+ 

dt 
dJ. 

dt 



Tb+ + T. 



p+ ■ 



Tb- + Tp 



(109) 
(110) 



By adding the expressions in Eqs. ( |109b and dl 101 ). we derive an 
equation for the time evolution of the global rotation Q^: 



dn, 

dt 



he Jo 



sin^e cos^ddO 



(111) 



By substracting Eq. ( II 10b from Eq. ( 11091 ) we derive after some 
algebra: 



ddO., 
dt 



(BlR^ r*.. ^ , ^ 
- — sin'*6' cos^^t/e 
he Jo ) 



(112) 



A characteristic damping time r appears to be defined by: 



T - 



ch 



= 3 X 10^ /45B;^U''^rosec- 



(113) 



When the magnetosphere is completely open, Eqs. dl 1 ll )-( fTT2b 
reduce to: 



dt 
dt 



2 

157 

21, 



Q.. 



■ 



(114) 
(115) 



35 /2T 

From Eqs. (ll04l )-( fT05] l. we find that the damping times for 
and in the open regime are, respectively: 



"^brake 



Todd - 



15t 



7.5 T, 

(■kh 



15^ sin'^ecos^e tan (1 cose) dO 



2n 



£'^sin^0cos^ed0 



= 7.: 



(116) 
(117) 



When the fast-spinning aligned rotator experiences episodes of 
magnetospheric opening, its evolution consists of a succession 
of open (or high) states and closed, classical (low) states. During 
open periods, Eqs. ( II 14b -( fTT5l ) apply and the energy output, of 
the order indicated by Eq. ( |92l ), is considerable. During these 
periods, the open field should occasionally reconnect, attempt- 
ing to return to a closed structure, but, once reformed, the latter 
is again blown open after the very short time needed to build 
a twist again of approximately half a turn. Large irregular vari- 
ability is then expected during these open periods, down to the 
millisecond timescale, which is the time to cross through a light- 
cylinder size, expected to be representative of the equatorial cur- 
rent sheet, at the speed of light. We note that the closed episodes 
are initially short in duration when Oodd is of the order of 10""* 
as supected. When a total time of order of T„dd has been spent 
in the open state, the oscillation has damped to an amplititude 
insufficient to open the magnetosphere, and the average rotation 
has been substantially reduced. Neglecting the weak damping 
experienced during the closed episodes, the spin rate Q*(f) and 
the largest amplitude of the differential rotation (JQ,(f) vary as 



(118) 
(119) 



where t is the cumulated time spent in the open state. The emit- 
ted power scales as QJ, and declines in a time of approximately 
15t/4. The GRB would disappear from view in about this time, 
which, for a polar field of 3 x lO''* gauss, is about 20 minutes. 
The openings completely cease when the maximum twist falls 
below half a turn, that is, from Eq.([86]l, when 5Q,/Qr ~ tt . This 
happens after a time ts,op such that 



' stop 



Todd In 



nQ.T 



Todd In 



y/aodd Pt 
ttPm 



(120) 



The mean rotation rate 0.tsiop at time ts,op is given by Eq. dl 18b . 
We define ^ = ToddlTtrake, which, from Eqs. dl 16b -( fTT7] i. is 
nearly unity. Relating 5Q,o to aodd as in Eq. dSTl i. we obtain, 
for ^ : 



Stop 



P* 



stop 



30 ms 



Ctodd 



y/aodd Pt Numerically: 

1/2/ Pj 



(121) 



For aodd = 3 X lO""* and Pt = 0.8 s (i.e. half the period ^ for 
Bp - 3 10'^ gauss) Ptstop ~ 4 ms. After spending a time ts,op in 
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the open state, the high episodes cease and the object becomes a 
pulsar with a period of the order of 5 milliseconds. As indicated 
above, the time of high activity is about 3.5 t (Eq. ( 1113b ). which 
is comparable to the observed timescale of long duration GRBs 
when the field of the compact star is somewhat higher than 10'** 
gauss. For high fields, of order lO'^ gauss, this timescale is about 
100 seconds. 



6. Conclusion 

It is natural to consider that a new-born quark star experiences 
differential rotation, causing its internal wound-up toroidal field 
to increase in strength to about 10'^ gauss. This motion then de- 
velops into a magnetic torsional oscillation, which could be the 
origin of long-duration y-ray bursts. We have indeed shown that 
an odd oscillation of small amplitude, which should be easily 
reached, is sufficient to open the star's magnetosphere. A similar 
effect would also result from other causes of north-south asym- 
metries. The rapid rotation then drives a relativistic wind from 
the entire stellar surface. When the star is a quark star, this wind 
is entirely leptonic. We have calculated the Poynting power re- 
leased and the timescale of this phenomenon, which meet the 
observational constraints if the polar field of the quark star is of 
the order of a few lO'** gauss and its initial rotational angular ve- 
locity is of the order of 300 Hz. Large amplitude variations in the 
light curve on timescales ranging from minutes to milliseconds 
is a natural outcome of this process. 
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Appendix A: IVIagnetosphere opening: an example 

We describe the asymptotic properties of the solutions to Eq. 
(l84l i in the context of a self-similar model and show that there 
is a limit twisting for closed solutions to exist. We define A to 
be the equatorial value of the flux function a (Eq. (f74ll). When 
a field line on the magnetic surface a is twisted, there is a rela- 
tion between its twist ifr(a), which is the difference in longitude 
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between its footpoints, and the poloidal current 1(a). The differ- 
ential equation of a field line is indeed: 



dr 



rd9 

"s7 



rsinOdip 



(A.l) 



Using Eqs. (f74l i and ( |79] |. we evaluate the change in longitude 
x}/ accumulated following a field line on the magnetic surface a 
from one of its footpoints Pi to the other P2: 



if/(a) - 1(a) 



f 

Jp, 



dtp 



r^sin^l Bp 



(A.2) 



where dtp is the line element along the poloidal field line a. 
It is the twist i/r, not the poloidal current /, which is known. 
The relation Eq. ( IA.2b does not infer 1(a) directly when ilj(a) 
is known, because the line of constant a on which the integra- 
tion in Eq. ( IA.2I ) is to be carried is unknown before the prob- 
lem expressed in Eq. ( l84b has been solved. To determine un- 
der which conditions a dipolar field line closing at a few stellar 
radii would be inflated sufficiently by the magnetospheric cur- 
rent to reach the light-cylin der, we attemp t to identify separable 
solutions to Eq.dgU i (Lvnde n-Bell & Boilv 1 994: Wolf son 1995; 
iBardou & Hevvaertslll996() of the form: 



a(r, 9)^ A 



(7) ^^'^ 



(A.3) 



IWolfsonl(ll995h numerically studied the solutions of Eq. (|84] l un- 
der the ansatz (|A.3l l. We show here that the solutions must open 
when the twist reaches a finite value. The largest value of a being 
the total star flux A, g(6l) < 1, andg(7r/2) = 1 . Similarly ^(0) = 0, 
since there is no flux through a circle of zero radius centred on 
the polar axis. The apex of field line a is at a distance D(a), such 
that D(a) - R(A/ay^'' : the more inflated the magnetosphere, 
the smaller the parameter p. We will then deal with the small 
p limit. A solution of the form of Eq. (IA.3I I cannot match any 
given flux distribution on the star, nor any given twist i//(a). The 
constraint of Eq. dAT2] i can only be satisfied in an average sense. 
Using Eq. (IA.3b in Eq. ( l84l i the following equation is obtained: 



pip- 



1) sine d 
— + 



1 



A^'P I(a)r(a) 



,1+2/;; 



(A.4) 



The left-hand side is a function of 6, while the right-hand side is 
a function of a. Both then equal a common constant, -K. This 
means that for the ansatz of Eq. iA.3i to be satisfied when the 
similarity exponent is p, the current function 1(a) must be: 



1(a) _ Kp fli+'^P 
~ ~ \pTT A^'PR 



(A.5) 



The dipolar angular function g - sin 6 is recovered for p - I 
and / = A" = 0. For non-vanishing K, the constraint of Eq. jA.H 
becomes, for solutions of the form Eq. ( |A.3l i: 



m 

c 



\p(a) 



A^iPR £gi/p(9)de 



(A.6) 



The limits on the integral at the denominator reflects the fact 
that all field lines span the interval [Q,n] in 6 (Eq. (IA.3I )). 
Equation ( |A.6t is consistent with Eq. (|A.5l l only when tfr(a) is 
independent of a, which corresponds to the peculiar twist profile 
in which one hemisphere rotates like a solid body and the other 
in an opposite sense. We define < > to be this twist. We may 



think of < (/r > as being some average on one hemisphere of a 
more realistic twist profile. K is related to < tfr > by: 



Kp 



P <^> 



p + i ^^g'ip(e)d9 



(A.l) 



Having chosen a value of < ip > this relation infers K for a given 
p. Given this relation, the value of p itself results from the need 
for the solution of the angular function g(ff) to satisfy the three 
requirements g(nl2) - 1, g(Q) - 0, and g'(nl2) - 0, the last 
one resulting from the symmetry of magnetic surfaces with re- 
spect to the equator. A second-order differential equation accept- 
ing only two boundary conditions, the extra condition eventually 
determines the value of p. To express this condition explicitly, 
Eq. ( IA.4b for the angular function g must be solved in the limit 
of small p. In terms of the variable x - cos 6, Eq. ( IA.4l i can be 
written: 



(\-x^)'^+p(p+\)g + Kg'^^IP 



(A.8) 



For small p the second term of Eq. ( IA.8I 1 is negligible. The expo- 
nent (1 +21 p) being very large, the third term of dA.St essentially 
vanishes wherever g < 1 . It remains non-negligible only in the 
vicinity of the equator (x - 0), where g reaches unity. The so- 
lution g(x) is then almost a linear function for all x, except in a 
small region about x = 0. This allows us to simplify Eq. (IA.8I) in 
the small p limit as: 



dx^ 



+ Kg 



1+2/ p 







(A.9) 



This equation has a first integral. The condition that g'(0) = at 
X - 0, where g must equal unity, can be satisfied by an appropri- 
ate choice of the integration constant, giving: 



_Kp 

P + 



(A. 10) 



Wherever g is sufficiently less than unity, g k Kp. However, 
we know that in these regions the modulus of the slope should 
be unity, because g has already been recognized to be a linear 
function varying from g = Q at x = 1 to very nearly ^ = 1 at 
X - Q. The relation between K and p is then, in the small p Umit: 



Kp^l ■ 



(A. 11) 



We should now establish the relation between < iff > and p re- 
sulting from Eq. (IA.7I ). To calculate the angular integral at the 
denominator, we must solve Eq. ( lA.lOt when p « 1 . Taking Eq. 
( lA.llt into account, changing the unknown function g for h such 
that g = (I - ph) and making use of the fact that the limit for p 
approaching of (1 + pyY^P is expCy), it can be shown that the 
solution of Eq. dA.lOt is in this hmit: 



g - 1 - /:> In cosh 



'(;)) 



(A. 12) 



Since p is small, g'^'^ » 1/ cosh(x//?). The integral that appears 
in Eq. ( IA.7l i can then be calculated, resulting in: 

(A. 13) 



lim < d/ >- n 



Thus, the exponent p approaches zero as the twist approaches n 
and the magnetosphere swells boundlessly in this limit. For our 
purpose, it is sufficient that a field line extends farther than the 
light-cylinder for it to open. We can then safely adopt the limit 
of a twist of a half a turn in causing an almost complete opening. 



